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We present detailed analysis of topological properties of the valence bond solid (VBS) states doped with 
fermionic holes. As concrete examples, we consider the supersymmetric extension of the SU(2)- and the SO(5) 
VBS states, dubbed UOSp(l|2) and U0Sp(l|4) supersymmetric VBS states, respectively. Specifically, we in- 
vestigate the string-order parameters and the entanglement spectra of these states to find that, even when the 
parent states (bosonic VBS states) do not support the string order, they recover it when holes are doped and the 
fermionic sector appears in the entanglement spectrum. This peculiar properties are discussed in the light of the 
symmetry-protected topological order. To this end, we characterize a few typical classes of symmetry-protected 
topological orders in terms of supermatrix-product states (SMPS). From this, we see that the topological order 
in the bulk manifests itself in the transformation properties of the SMPS in question and thereby affects the 
structure of the entanglement spectrum. Then, we explicitly relate the existence of the string order and the 
structure of the entanglement spectrum to explain the recovery and the stabilization of the string order in the 
supersymmetric systems. 
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I. INTRODUCTION 

The valence-bond solid (VBS) states has been originally 
introduced by Affleck, Kennedy, Lieb and Tasaki' to build 
explicit model ground states which realize the properties of 
the generic integer-spin antifeiTomagnetic spin chains conjec- 
tured by Haldane.^ Quite unexpectedly, on top of the proper- 
ties already anticipated from other analyses (e.g. quantum- 
disordered ground state with short-range spin correlations, 
gapped triplet spin excitations, etc.), these states exhibit many 
striking features such as the emergent boundary excitations 
(edge states jS and the existence of hidden string order.SEl 
In the case of spin-1 systems, it has been arguecP that the 
hidden topological (string) order is a consequence of the 
Z2 xZ2-symmetry breaking occurring in the system after ap- 
plying the non-local unitary transformation. The idea of non- 
local hidden order and edge states have been to some extent 
generalized^ '° to other values of integer-spin-S" although the 
hidden Z2 XZ2 -symmetry is never broken'^ in the case of even- 
S. Through these studies, it has been recognized that there are 
some differences^ ^ in the ground-state properties according to 
the parity of S. Nevertheless, by analogy with the quantum- 
Hall systems", the ground state of generic integer-spin anti- 
ferromagnetic chains, including the original VBS state and its 
higher-spin generalizations'^, characterized by certain kinds 
of non-local congelations and emergent edge states have been 
called 'topological' in a rough sense. 

Recent development in quantum-information-theoretic ap- 
proaches to quantum many-body problems enables us to ex- 
tract information on the bulk topological order from the en- 
tanglement properties of the ground-state wave functiorfSE] 
The topological states in one-dimensional spin systems have 
been reconsiderecS^CHfrom the modern point of view and the 
precise meaning of the topological Haldane phase has been 
clarified. In these studies, the string order parameters and the 
edge states, which in general are not robust against small per- 
turbations, are replaced by more robust objects (i.e. the struc- 
ture of the entanglement spectrum or the structure of tensor- 



network). In particular, it has been shown in Ref. TTthat the 
existence of (at least one of) the discrete symmetries -time- 
reversal, link-inversion and Z2 x Z2 symmetry-divides all 
states of matter in ID into two categories-topologically-non- 
trivial ones and the rest; generic odd-integer-5 spin chains 
belong to the former while even-S* chains to the latter. The 
hallmark of the topological phase protected by the above dis- 
crete symmetries is that all entanglement levels are even-fold 
degenerate. In this formulation, the difference between odd-S" 
and even-S* is naturally understood in terms of the entangle- 
ment structure. 

In this paper, we present an exhaustive discussion about 
the effects of coexisting bosonic- and fermionic degrees of 
freedom on (symmetry-protected) topological phases in ID. 
Clearly, this kind of questions is motivated in part by hole 
doping in the Haldane-gap systems.'^ ^" In order to incor- 
porate the coexisting bosons and fermions, for mathemati- 
cal convenience, we use supersymmetry (SUSY) which re- 
lates bosons caiTying integer spins and fermions with half- 
odd-integer spins. Several 'topological phases' with SUSY 
have been found so far in, e.g. quantum-Hall systems^', VBS 
state and ultra-cold atom systems^-. However, the pre- 
cise characterization of these SUSY topological phases has 
not been obtained so far and it would be quite useful to inves- 
tigate symmetry-protected topological order in model SUSY 
systems from the entanglement point of view. 

As the model SUSY states, we consider a class of super- 
symmetric VBS (SVBS) states defined by the Schwinger op- 
erator consisting of 2K bosons which represent the bosonic 
degrees of freedom at each site (e.g. localized integer spins) 
and N fermions which correspond to doped fermionic holes 
(with K and being integers). This class is interesting since 
it includes the SVBS states investigated in Refs. .22, and l24l 
as well as the SUSY-extension of the SO(5) VBS state and 
the Sp(A^) VBS state introduced respectively in Ref. 'SS'and 
Ref. |26l The (S)VBS states are rare examples where we can 
study non-trivial topological properties even in ID and most 
of calculations can be done without relying on any approxima- 
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tion. Taking advantage of such properties of the SVBS states, 
we uncover the roles of SUSY in topological phases in ID. 

The generalized hidden string order in the SVBS stateJ^ 
has been investigated already in the previous work^"* by the au- 
thors. In contrast to what is known for the bosonic counterpart 
(the spin-S* VBS statJ^, the symptom of the non-trivial topo- 
logical order has been observed in the analysis of the string 
order even for the even-integer superspin. To be more precise, 
even when the string order vanishes, it revives upon the hole 
doping; this might suggest the existence of topological order 
in the SVBS states regardless of the parity of bulk superspin 
S. In order for the better understanding of this phenomenon, 
we first characterize symmetry-protected topological orders 
in SUSY systems in the language of entanglement. To this 
end, we use the supermatrix-product-state (SMPS) formalism 
to generalize the arguments of Ref. JLZand derive the relation 
between topological order in the bulk and the entanglement 
structure. The SMPS formalism further enables us to obtain 
the explicit relation between the entanglement spectrum and 
the string order parameters, and thereby to clarify why the 
hidden string order revives after doping. 

As has been emphasized in the previous workP^ in spite of 
its name, the SMPS formalism does not assume any particular 
form of SUSY. In fact, we do not need even postulate exact 
SUSY and the only prerequisite is that the local Hilbert space 
is made up of the bosonic part and the fermionic one. In view 
of the ability of (S)MPS in approximating any gapped states in 
ID with arbitrary precisiorl22l2SI our results are applicable to 
a wider class of ID systems with some kind relation between 
bosons and fermions. 



The organization of the present paper is as follows. In sec- 
tionjll] we introduce a class of UOSp(iV|2ii')-invariant SVBS 
states (2K being the number of boson species and N for 
fermions) with arbitrary superspins using the Schwinger op- 
erator. We then construct the explicit SMPS representation 
for {N, K) = (1, 1) (UOSp(l|2)) and (1, 2) (UOSp(l|4)) and 
summarize several important properties of these states. As 
the first step toward the investigation of topological order, we 
explicitly evaluate the string order parameters in the above 
two types of SVBS states for different values of superspins 



in Section III There we find that the revival of the string or- 
der already observed for UOSp(l |2) in Ref.|24]occurs in other 
SUSY cases as well. In section|lV] the entanglement spectrum 
of these SVBS states (in the limit of infinite-size systems) is 
derived and typical features of the spectrum is discussed. In 
order to understand the results obtained in the previous sec- 
tion and characterize symmetry-protected topological order in 
ID SUSY systems, we generalize the argument of Ref.lTTlto 
SUSY systems in section|V]and relate the structure of the en- 
tanglement spectrum and the bulk topological order Finally, 
the relationship between the degeneracy of the entanglement 
spectrum and non-vanishing string order parameters is clari- 
fied in section VI by using the (S)MPS formalism. Section 



II. SVBS STATES AND SMPS FORMALISM 

In this section, we briefly describe how the standard MPS 
formalism is generalized to the cases with SUSY. Let us be- 
gin with constructing the MPS of the spin-Af {M: integer) 
SU(2) valence-bond solid (VBS) state'^ starting from its rep- 
resentation in terms of the SU(2) Schwinger operators (f) — 



|VBS)(^^)=n(&f^\/-fef&]+J)'^|vac 
i 

= n('^^-*'^2<^:'+i)*'|vac), 



(1) 



where the metric (or, the SU(2) charge conjugation matrix) 



ia2 



1 

-1 



(2) 



has been used to form a maximally-entangled (singlet) pair 
between the site j and j + I- Therefore, by construction, the 
VBS state is SU(2) invariant and represents a spin-isotropic 
state. 



A. General Idea 

The standard construction of the VBS-type of stateJ^starts 
by preparing two auxiliary degrees of freedom on each site of 
the lattice. Then, the (bosonic) VBS state is constructed first 
by creating singlets between pairs of those auxiliary objects 
on adjacent sites and then by projecting the tensor-product of 
the two auxiliary objects on each site onto the desired physical 
Hilbert space. 

The SVBS states are introduced by including the states 
with one- or more fermionic holes into the above bosonic 
Hilbert space. Mathematically, we replace the usual Lie- 
group symmetry (e.g. SU(2)) with that of the super Lie group 
UOSp(A^|2iir) corresponding to 2K bosonic degrees of free- 
dom and N fermionic ones (for a review of super Lie groups, 
see, for instance, Ref. 30, and for UOSp(iV|2/-i:), Ref. 31). 
SpecificaUy, the SVBS states with UOSp(iV|2i4:)-symmetry 
is defined as 

|SVBS(7V|2X))(^'^) = n(^J^A^|2/^V',)'^|vac) , (3) 
where ijj stands for the UOSp(Af |2iir) Schwinger operator 

^^{b'\b'\---y'<\f'\---,r'y. (4) 

The 2K bosons 6'^^ (o- = 1, 2, • • • , 2K) and the N fermions 
(/i = 1,2,- •• ,N) satisfy the commutation relations 

= j.r^ {f'',r^} = [b'',r] = [6^,/'^^] = o. 

The matrix 7?.jv|2if signifies the \JOSp(N\2K) invariant ma- 
trix: 



VII is devoted to summary and discussions. 



7^ 



N\2K 



J'. 



2K 





-In 



(5) 
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where the USp(2i4r)-invariant 2Kx2K antisymmetric matrix 
J2K is defined using the PauH matrix (T2 as: 



J: 



/«CT2 \ 



2K 



\ 



(6) 



1(72 



and In denotes the A^-dimensional identity matrix. By 
using the above equations, it is straightforward to show 
that the product of spinors '4'iT^N\2K'^j is singlet under 
UOSp(iV|2/-i:). 

As the number of fermion species TV corresponds to that of 
the SUSY in the system, hereafter we call the SVBS states 
defined by ^ and @ the UOSp(Ar|2i^) SVBS states. In the 
present paper, we give the detailed discussions for the two 
N ^ 1 cases, specifically iK,N) = (1,1) and {K,N) = 
(2, 1), in which the following isomorphisms between the 
orthogonal groups and the unitary symplectic groups hold: 
SO(3) ~ USp(2)/Z2 (K = 1), SO(5) ~ UOSp(4)/Z2 
(K = 2). For UOSp(iV|2) (K = 1), the metric matrix is 
given by 



■77 _ {i'^2 

l<-N\2 - 1 

and for UOSp(7V|4) {K ^ 2), by 





'In 



n 



N\i 



'ia2 
ia2 
-In, 



(7) 



(8) 



The particle number at each site is related to the superspin 
S via 



25 



2K+N 
^ /,"t,/,Q 



2K 



N 



a—1 (T—l 



b^h" + f^f" = zM, (9) 
^=1 



where z is the lattice-coordination number (z = 2 in one di- 
mension). Throughout this paper, we reserve the symbol S for 
superspin and use S for the bosonic spin. Since f^^f^ 
takes either or 1, the possible values of SU(2) spin, which is 
equal to the half of the number of bosons at each site, are: 



1 



2K 



2 ^ 

= \zM, \zM - \zM - 1, ■ 
2 ' 2 2' 2 



-zM --N. 
2 2 



(10) 



(If N > zM, it is implied that the above sequence terminates 
at 5 = 0). One may find that the inclusion of SUSY intro- 
duces, as well as the states with the spin magnitude zM/2 
which exist already in the SU(2) case, those with spin smaller 
by 1/2. In what follows, we consider the one-dimensional 
cases (i.e. z = 2) unless otherwise stated. 

For the ID chain [z — 2), the above sequence reads 



S = M, M 



M - 1, 



M 



1 



(11) 



and correspondingly the emergent edge spin takes the follow- 
ing values 



s = -M, -M- ^, 
2 ' 2 2' 2 



1, 



2 2 



(12) 



(again, if > M, the above sequence is understood as to 
stop at s ~ 0.) The dimension of the physical Hilbert space 
at each site constructed in this way is given by the sum of the 
one of each bosonic Hilbert space with a fixed boson number 
(25 - n): 



N 



dsiN\2K) = 



ri=0 



f2K + 2S-n-l 
y 2K -I 



(13) 



It should be noted here that the Sch winger-operator construc- 
tion presented here does not cover all the possible VBS- 
type states with UOSp(A^|2_ftr)-s^mmetry. In fact, there is 
an important class of VBS states-^^ which is a SUSY gener- 
alization of a series of SO(2ri + l)-invariant and USp(2if )- 
invariant states considered respectively in Refs. 33 and 34 and 
in Ref. 26, However, most of the conclusions obtained here 
hold for those models as well. 

The UOSp(7V|2A') SVBS state ^ may be rewritten as 

|SVBS(A^|2if))(^^) = n(^'^^|2/^V'.+i)'''|vac) 

i 

^X{mR%l^^^+l)Wc), (14) 

i 

where is a graded fully symmetric representation of 
UOSp(A^|2A:) of the order M and R^n\2K metric for this 
representation.^' Another equivalent form (a matrix-product 
form)^"^ may be useful for practical purposes: 



|SVBS(Ar|2i^))(^'^) = • • -^i, 

where the matrix Ai is defined as: 

A, ^ <g^*,*',|vac), . 

B. UOSp(l|2) SVBS states 



(15) 



(16) 



Let us begin with the simplest casd^(Ar,A:) = (1,1). 
The graded Schwinger operator is given by 



(17) 



and the corresponding SVBS state, which we call the 
UOSp(l |2) SVBS state (precisely, this is the one dubbed type- 
I in Ref. 24), is given by: 

|SVBS(1|2))(*^) =n(46l+i-6l«li-^///i-i)''|vac) , 

(18) 

where we have added the fermion doping parameters r by 
hand. However, such a parameter may be absorbed in the re- 
definition of the normalization of fermions (/^ 1— > P l\Jr, 
f I-?- v^/) and the SVBS states possess the SUSY even for 
finite values of the parameter r. 
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1. S=l 

Let us consider the superspin 5 = 1 case. Since S is related 
to the number M of SUSY valence bonds through (|9|, the case 
M = 1 of eq.([r8|) corresponds to 5 = 1. 

The SVBS state on a finite open chain is specified its edge 
states, a and /3, respectively on the site 1 and L: 

L-l 

|svBS(i|2))W = (7^l|2Vl)" n(^^^i|2^^+i) ^^\^^^)' 

1=1 

(19) 

where i/'} = (aj, and the UOSp(l|2) metric 7^l|2 

is defined in Jvl. The state ISVBS-I)!^^"^'' can be expressed 



as a product of the matrices defined on a each site: 



|SVBS(1|2))« = (4^)4^^ 



(20) 



where 4^'' is given by 

4^) = 7^f V,^>ac), 

|0), ^/2| - 1), 
-V2\l), -\0), -^^|l/2), 
,^V^\l/2), -V^\-l/2), 

a=-l,Oa <T=-l/2,l/2 

with 

A{-1) 

\0 


A(l/2) =1 -Vr I , A(-l/2) 

V?^ 



(21) 



0\ 







V2 




-1 


Oyi 







V2 0\ 













oJ 







The five basis states corresponding to the 5 = 1 irreducible 
representation (denoted by 5) are given by 

|1) = i=at'|vac), |0)=at&t|vac), | - 1) = -^6t'|vac), 

|l/2) = aV^lvac), | - 1/2) = b^p\wac) , (23) 

where |vac) is the vacuum of both the boson and the fermion: 
a|vac) = 6|vac) = /|vac) = 0. The first three states cor- 
responds to the spin-1 (3) representation of SU(2), and the 

second two states constitute 2 with spin-1/2. 

The parent Hamiltonian of the state ([19]) is constructecP^I^ 
in such a way that the local Hamiltonian /ij j+i acting on the 
bond (j, j + 1) annihilates all the nine states appearing in the 
product AjAj+i- Therefore, the ground state on a finite open 
chain is nine-fold degenerate with respect to the matrix in- 
dices. Since the i/jj and V'} represent the two auxiliary degrees 




of freedom at the site j, the above nine-fold degeneracy re- 
flects the existence of the three edge degrees of freedom on 
both edges of an open chain: 

It)) =a^ I vac), i;)) =6t|vac), |0)) = /t|vac) . (24) 

As the doping parameter r is changed, the state ( [T9] l in- 
terpolates between the two well-known states: at r — >^ 0, 
|SVBS(l|2))(i) is reduced to the original VBS stat^|VBS) 

|SVBS(l|2))(i) ^ |VBS)(i) =n(«l&I+i -fel«m)|vac), 

(25) 

while, at r — > oo, |SVBS-I) is reduced to the Majumdar- 
Ghosh (MG) dimer state^^ |MG) 



|SVBS(l|2))W^n//|MG), 



(26) 



where 



|MG) = ( n - \{)i^^A+i - &lali)|vac). (27) 



z:even i;odd 



In the discussion of the entanglement spectra (section|IV|i, we 
will see in the two limits, the entanglement entropy nicely in- 
terpolates between that of the VBS state and the MG state. 



2. Higher-S 

It is easy to generalize the above strategy to the cases with 
general superspin-5. In Ref. |24l the expression of the A- 
matrix for superspin-5 type-I SVBS state is given as: 

= ^.^(a], 6], /j) J-'^(a], 6], /;)|vac), , (28) 

where the 5-th order polynomials J'^ and are defined in 
eqs.(C3a) and (C3b) of Ref.^ The above expression may be 
readily rewritten into the standard form ( [T6| ): 



4tnj) = <2^*,*;|vac),, 

where 

{^,)a^T^{a] , b] , /j) (l<a<25 + l) 



(29a) 



_ (-!)' 



3fc,(5+2)-a 



(1 < a,5 < 5 + 1) 



(-lf-i-i)S,^f^3S+3)-a (5 + 2 < a, 6 < 25 - 



-1), 
(29b) 



C. UOSp(l 4) SVBS states 

Now we proceed to the case (A'', K) = (1, 2) (one 
fermion species and four bosonic). For UOSp(l |4), the graded 
Schwinger operator is given as: 



,p = {b'\b'\b^\b^\V-rfy 



(30) 



5 



These five operators correspond to the five-dimensional rep- 
resentation (5) of UOSp(l|4); the first four ib^\ b^'' ) 



respectively create the four bosonic states 



|1) 
|3) 



1 1 

2' 2 
1 1 

2'^2 



, |2) = 
|4> = 



1 1 

2'"2 
1 1 
^2' 2 



(31) 



which are already contained in the spinor representation of 
SO(5) and the last one creates the fermionic state 1 5) — \f). 
We prepare z copies of 5s to construct the physical Hilbert 
space at each site of the lattice with the coordination number 
z and, according to which representation is chosen from the 
tensor product of z 5s, we can obtain several different types 
of MPSs. For instance, since a pair of 5s is decomposed as 



5(Ki5 - 10 10 14, 



(32) 



two different SVBS states (10 and 14) are obtained in one 
dimension (z = 2). 

Following the general method described in section II A one 
can construct the following UOSp(l|4) SVBS state; 



|SVBS(1|4))(*^) = n(^'7ei|4V,)*^|vac) 



{id) 



ltl,2t ,2t,lt , ,3t,4t ,4t,3t 



Kb, 



K^ 



K^ 



(33) 



where the summation is taken over the nearest-neighbor pairs 
(i, j) and r denotes a real parameter varying from to cx). 
The state has the same structure as the UOSp(l |2) SVBS state 
except for the metric Ti-m defined in (|5]) or ([8]). The superspin 
S in this state is given as 



2S^Y.^lhl+f}h = zM 



(34) 



The dimension of the local physical Hilbert space (i.e. the size 
of the representation S) ^ reads for {N, K) ^ {1,2): 



(45 + 3)(25 + l)(5 + l) 



(35) 



In the following, we consider the one-dimensional case (z = 
2) with M — 1{— S) where the SO(5) spin magnitude takes 
the following two values: 



1 4 . 



(36) 



and di(l|4) = 14. 



On a finite one-dimensional chain, the UOSp(l|4) SVBS 
state ((33]l may be written as 



|SVBS(T)), 



L-l 



= {^i|4^i}„, n(^}^i|4^J+i) {^Uc, Kac) (37) 



(T) .(T) 



l(T) 



where 7^i|4 is given by (jSjl with = 1. The matrix A is 
defined by 



/ |1,2> 
-%/2|l,l) 
|1,4) 
-|1,3) 



V2|2,2) 
-|1,2> 
|2,4) 



|2,3) 
-|1,3> 
|3,4) 



-|2,3) -V2|3,3) 



|2,4) 
-|1,4) 
V2|4,4) 
-|3,4) 



V-^|l,/) -^A^|2,/) -V^|3,/) -^/f|4,/) 



%^|2,/) \ 

v^|4,/) 
-V^|3,/) 
/ 



4 

E 

a<T=l 



Af(a,r)|a,r) + ^<>(a)|a,/), 



where the D ^ 14 basis states are given in terms of the graded 
Schwinger operators in (l30| as (a, t — 1,2, 3, 4): 



1 



\a, a) 



^/2 



<Tt^2| 



(38) 



|ct,t) = 6'"^6^Vac> (o- < r) , 
h/>^ 6- Vulvae) . 

The expressions of the 14 matrices A{(t, t) and A{(j) are 



given in appendix A 1 



Since the Schwinger operators are used, it is obvious that 
the physical Hilbert space thus constructed is the 5 = 1 (i.e. 
14) fully symmetric representation in the tensor-product de- 
composition ([32]i: 



Jfully-sym. 



14 



S0(5) 



> 10 4, 



(39) 



where >' denotes the decomposition into the SO(5) ir- 
reducible representations. As in the case of UOSp(l|2) 
{{N,K) — (1, 1)), the physical Hilbert space contains two 
irreducible representations of SO(5): the spinor- (4) and the 
adjoint (10) representations. Since all the 14 basis correspond 
to the components of the rank-2 symmetric tensor made of the 
two constituent spinors (5), we call the MPS thus constructed 
tensor-type. 

A remark is in order here about other possible MPSs. In 
fact, as has been mentioned before, another important MPS is 
obtained-'^ if we use the 10-dimensional anti-symmetric repre- 
sentation (vector representation; hence the MPS may be called 
'vector- type'), in stead of the 14-dimensional one 



S0(5)_ 



(5 «) 5)a„ti-sym. = 10 — ^ 5 4 1 



(40) 



The MPS obtained in this way is a direct generalization of 
the SO(5)-invariant MPS considered in Refs.[33]and|34l The 
details of this class of MPS will be reported elsewhera^^. 
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1. Limiting Cases 



A. UOSp(l|2) SVBS states 



Now let us consider the two important limiting cases r — > 
and r ^- oo. In the limit r = 0, the UOSp(l|4) SVBS states 
( |33| l or ( (37] i reduce to the following VBS states 

|VBS> = n &f - &f + &f fef - V)''|vac), 

(41) 

dubbed bosonic S0(5) VBS state in Ref.'ZS^ 

In the other limit r —i' oo, the dominant part of yl*^* reads 
(after dropping factors proportional to ^/r) 



/ 






















































l 


-ll>, 


-|2), 


-|3), 


-|4), 



|2>, \ 
-|1>/ 
|4>, 

-|3>, , 
/ 



(42) 



Then, the two-site MPS + takes the following 

block-diagonal form 



± 



s2.2 0-,j + i); ' 

(43) 

where |1) , . . . , |4) are defined in eq.((3T| and the (2, 2)-block 
is the SO(5)-singlet made up of two spinors: 

S2,2(j,J+l) - |l),|2),+i-|2),|l),+i + |3>,|4),+i-|4),|3),+i 

(44) 

When the 4x4 matrix Bi^i{j, j + 1) is multiplied by Bi^i{j + 
2,j + 3) from the right, a new SO(5)-singlet is inserted at 
the bond (j-|-l,j+2). Therefore, one sees that the string 
of represents an SO(5)-generalization of the Majumdar- 
Ghosh valence-bond crystal. The vector-type UOSp(l|4) 
SVBS state mentioned above shares the same property.-^ 



state-A • 



maximally-entangled 



1 



state-B 



FIG. 1. (color online) The r oo limit of 5 = 1 SVBS state. Filled 
circles denote the bosonic qubits (S = 1/2 spins for U0Sp(l|2) 
and 4-dimensional SO(5) spinors for U0Sp(l|4)). On a chain with 
even number of sites, the MPS is block diagonal with the (l,l)-block 
Bi^i and the (2,2)-block 62,2 corresponding to state-A and B, re- 
spectively. 



III. STRING ORDER 

One of the striking features of these VBS states is the exis- 
tence of non-local order called string order. In the usual spin 
systems, it is known^ that the string order is a manifestation 
of the spontaneous Z2XZ2 symmetry breaking in the ground 
state. 



In the case of the usual (pure) spin systems, the string or- 
der parameters are defined by the infinite-distance hmit of the 
string correlation functions^ 



OJtring = lim ( exp 



j+n-1 
^TT ^ SI 

k=j 

j+n 



S 



j+n 



S: 



j+n 



(45a) 



(45b) 



It is straightforward to generalize the string order param- 
eters to the case with SUSY by replacing the spin operators 
5° to their (45+1) -dimensional expressions. For superspin 
5 = 1, it is given bj^^ (OstrinR = O^aino by SU(2)-symmetry): 



stnng 



(r) 



4 {r^ + UP + 18 + 2 (7-2 + 3) V8r2 ^ 9} 
(8r2 + 9) (V8r2 + 9 + 3)^ 



(46) 

In the limit r ^ 0, the above string expression reproduces 
the well-known value'^4/9 (perfect string correlation). In the 
opposite limit r 00, the string order parameter O^.^^^ ap- 
proaches to a finite value 1/16, which implies that the string 
order survives in the r 00 limit. This agrees with the fact 
that the spin-1 Haldane state is adiabatically connected to the 
spin-1/2 dimer state.E^ 

One can readily generalize the above results to the higher- 
S cases, 2^ which are SUSY-analogues of the higher-spin 
(bosonic) VBS state introduced in Ref. 12i In the original 
spin-5 VBS states (r — 0), the string order parameters have 
been investigated^ ^ and it has been concluded that they vanish 
for even integer S. In contrast, for finite values of the doping 
parameter r, the string order parameters revive^'* due to the 
existence of SUSY (see Fig.[2]i. This interesting behavior will 
be discussed in section [Vl] in the light of symmetry-protected 
topological order. 




l + r2 



FIG. 2. (color online) The string order parameter O^i^g for several 
values of superspin S plotted as a function of r [Ref. 24 1. Note that 
O^iagir~0) = for even-5 corresponding to the vanishing of string 
order parameter for even-S'. 
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B. UOSp(l|4) SVBS states 

In Ref. [34] it has been pointed out that the idea of 
hidden-symmetry breaking^ and the associated string order 
parameters'* can be generaUzed to a class of models with 
higher symmetry SO(2n+l) by using the 2" -dimensional 
spinor representation as the auxiliary Hilbert space. 

The four string order parameters for the SO(5) (n — 2) 
VBS state are definecPby analogy with their SU(2) cousin: 



o: 



ah 



strmg - lim ( exp 

* n/'oo\ ■' 



T ah 



(47) 



(L"^ = —L'"^ are the SO(5)-generators). The set of integers 
(a, b) (with a,b = 1, 2, 3, 4, 5) labels the ten generators and 
we may choose e.g. (o, 6) = (1, 2), (2, 5), (3, 4) and (4, 5). 

Since, by the SO(5) symmetry, the string order parame- 
ters are independent of the SO(5) indices a, b, we can assume 
(a, b) — (1, 2) without a loss of generality. In Ref. l34l it has 
been argued that the string order of the SO(5) VBS state is 
a consequence of the hidden (Z2xZ2)^ symmetry breaking. 
In the original SU(2) case, we pick up a pair {S^, S'^} and 
the two commuting Z2S are generated by e*^'^ and e"'^ , the 
former of which plays the role of the flipping operator of S^. 
In the case SO(5), we have two(=rank of SO(5)) such pairs 
(e.g. {L^^, L^*^} and {L^"*, L'*^}) and this is why the squai-e of 
Z2 XZ2 appears. Similarly, as we already know that the gener- 
alized string order exists^'' in the UOSp(l|2) SVBS state, we 
can expect finite string order in the case of UOSp(l |4) as well 
by considering two pairs of string order parameters. 

First we set r = and consider the SO(5) limit. By plot- 
ting the eigenvalues of local {L^^, L^'^) appearing in the string 
( |37j i of A^^\ one can easily see^^ that both L12 and L34 ex- 
hibit a kind of hidden antiferromagnetic order which is es- 
sentially the same as that observed^ in the 5 = 1 VBS state. 
In fact, the string order parameters ( |47] ) for (a, &) — (1,2) 
((a, b) = (3, 4)) removes the effects of the randomly inserted 
zeros in the L^^ (L^"^) configuration to pick up the hidden an- 
tiferromagnetic order 

The generalization of eq.(|47]l to the UOSp(l |4) SVBS state 
with arbitrary superspin-5 is straightforward; for 5 = 1, the 
bosonic generators L"'' are replaced by the 14-dimensional 
matrices (the explicit forms of them are not very important). 
The MPS formalism enables us to obtain the following result: 



Qab ^ 
^ string 



{4r2 + 3 (Vl6r2 + 25 + 5) } 
25) (Vl6r2 + 25 + 5)' 
(r -> 0) 




(48) 



In order to highlight qualitatively different behaviors with re- 
spect to the superspin S, we plot the result in Fig. [3] together 
with that of the superspin-2 case 



Qab 
^string 



49 7- 



V40r2 



49)' 



400 (40r2 + 49) 



(49) 



From this plot, one can clearly see that, for finite doping, both 
the 5=1 state and the 5 = 2 one are topological, while the 
latter is non-topological (i.e. non-Haldane) at r = (see also 
Fig.|2|l. The limiting value 1/16 is equal to the string order of 
the 5 = 1 UOSp(l|2) SVBS at r ^ 00. Similar results have 
been obtainecP^for the vector-type MPS mentioned in section 
HTCl 




1 + r2 



FIG. 3. (color online) The infinite-distance limit of the string cor- 
relation function for U0Sp(li4) states. At r = 0, the string order 
reproduces the known resulc23 9/25 = 0.36. Also plotted is the 
string order of the 5 = 2 (M = 2) state. As in the U0Sp(ll2) case, 
the string order vanishes at r = and revives after doping. 



IV. ENTANGLEMENT SPECTRA OF SVBS STATES 

In the pioneering paper, Li and Haldane'^ argued that the 
entanglement spectrum, which is obtained by taking logarithm 
of the Schmidt eigenvalues (or, the eigenvalues of the reduced 
density matrix) of the ground-state wave function, might be 
the fingerprint of the physical edge states that reflect the topo- 
logical order in the bulk. Since entanglement cut creates point 
boundaries in one dimension, we may expect that the discrete 
level structure of the entanglement spectrum reflects the bulk 
topological order. 

In order to carry out the explicit calculation of the Schmidt 
coefficients (or, entanglement spectrum), we adopt the SMPS 
formalism introduced in our previous paper.^'* One of the 
biggest merits of using the SMPS formalism is that the 
Schmidt decomposition, which is the essential step of the 
calculation, is almost done already when we write down the 
SMPS expression. Therefore, all we have to do is to rewrite 
the SMPS into the form of the Schmidt decomposition by us- 
ing the singular-value decomposition.^^ However, when the 
(S)MPSs with different edge states are asymptotically orthog- 
onal to each other in the infinite-size limit (this is the case in 
all (S)MPSs discussed below), the entanglement spectrum is 
most easily obtained from the (infinite-size) norms for differ- 
ent edge states: 



A. 



lim 

j-L-j.L/^Qo 



jAfj (q!l, a)ML-j{a, Qr) 



(50) 
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FIG. 4. (color online) The behavior of entanglement spectrum of 
the 5 = 1 U0Sp(l|2) SVBS state (the inset is for the bosonic-pair 
VBS state). 'B' and 'F' denote bosonic- and fermionic part of the 
spectrum, respectively. 



where Mj is the squared norm of the MPS on a length-j sys- 
tem 

AG(a,/3) = |(A^2---A)c.,/3p . (51) 

A. UOSp(l|2) SVBS states 

1. 5 = 1 

By utilizing the SMPS, the Schmidt coefficients of the 
SVBS infinite chain, are readily derived as 



Af 



Ai^ = A 



1 



Af'^ = Aa^ 



4 4\/9 + 8r2 ' 
1 3 



2V9 + 8r2 ' 



(52a) 
(52b) 



which are shown in Fig. |4] and the corresponding entangle- 
ment entropy 



(53) 



is also depicted in Fig|5] From the entanglement spectra, we 
find that the bosonic and the fermionic sectors exhibit distinct 
behaviors. As mentioned in section |1V A| the SVBS chain 
interpolates the original VBS chain (r — 0) and the MG dimer 
chain (?■ oo). Then, we expect the entanglement entropy 
of SVBS chain also reduces that of VBS at r — 0, and that of 
MG at r OD. Indeed, in such two limits, the entanglement 
entropy gives those of the VBS chain and MG dimer chain: 



lim S^^(r) = log 2, 

r— fO 

3 

lim S^^{r) = - log 2. 

r— >-oo Z 

The states are maximally entangled when 

Ai^ = A2' = A3' = 1/3 (atr = 3), 



(54) 



FIG. 5. The behavior of the entanglement entropy of the 5 = 1 
UOSp(I |2) SVBS state. (The inset is for the bosonic-pair VBS state.) 



where the entanglement entropy takes the maximal value 
^Onax) ^ Xo'g'i. In contrast to the usual bosonic VBS 
states,™ 

the entanglement entropy S'ee of the SVBS states 
differs from what is expected from the dimension of the MPS 
matrices (i.e. bond dimension); they attain the maximal en- 
tanglement only at a particular value of the doping parame- 
ter r, which is different from the position of the maximal en- 
tanglement of the corresponding maximally-entangled pairs 
[for more details, see the supplementary online documents 
Ref-Ha. 

The 'level crossing point' (r = 3) between the bosonic 
spectrum and the fermionic one generally does not imply a 
quantum phase transition, in the sense that divergence of phys- 
ical quantities, e.g. spin-spin correlation length, does not oc- 
cur at the point. The (open) S — \ SVBS chain accommo- 
dates 5 = 1/2 superspins at the edges, i.e. the number of the 
edge degrees of freedom is 3 corresponding to |vac), fe^ |vac) 
and /^^|vac). Therefore, as has been found**" ^' in the usual 
bosonic VBS states, one sees that the entanglement entropy is 
bounded by the logarithm of the number of the edge degrees 
of freedom. However, here is one remarkable point; since 
the parameter r controls the contributions of the bosonic- and 
the fermionic degrees of freedom, one might expect that the 
entanglement is maximal at r = 1 where they appear with 
equal amplitudes (indeed, this is the case for a system of two 
5 = 1/2 superqubits [see Ref. 42|). Contrary to this naive 
expectation, the explicit calculation indicates that the maxi- 
mally entangled point is located at r = 3 due to many-body 
effect of SUSY. Note still in the bosonic many-body case, the 
entanglement is maximal at r = 1 (see the inset in Fig|5]l. 

To see a property peculiar to the SUSY states, let us intro- 
duce a "boson-pair VBS state": 

Ib-p.-VBS) = n(°]^]+i - - '^cjc]^, )|vac), (56) 



where cj denotes a bosonic creation operator that satisfy 



Ci,cJ] = 8ij and satisfies aja^ 



2. The new 



(55) state |b-p.-VBS) derived simply by replacing the fermionic 
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operator in the SVBS state ([TSj with bosonic one c' nei- 
ther has the inversion symmetry with respect to the center of a 
Unk (link-inversion) nor has the UOSp(l|2) symmetry. More 
importantly. The entanglement spectrum is plotted in the inset 
of Fig. |4] As in the 5 = 1 SVBS state, the boson-pair VBS 
chain has three Schmidt eigenvalues, two of which are doubly 
degenerate the other is non-degenerate. On the other hand, 
the entanglement entropy (see the inset of Fig. |5]) exhibits a 
different asymptotic behavior for r —i' oo since |b-p.-VBS) 
reduces, in the limit r — > oo, to the product state Yij c] |vac), 
while the SUSY version |SVBS(1|2)) still retains finite entan- 
glement due to SUSY. 
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FIG. 6. (color online) The entanglement spectrum and the entangle- 
ment entropy (inset) of the 5 = 2 UOSp(l|2) SVBS chain. 'B' and 
'F' denote bosonic- and fermionic part of the spectrum, respectively. 



Next, we proceed to the 5 = 2 SVBS chain. The bulk su- 
perspin is 5 = 2 which consists of SU(2) 5 = 2 and 5 = 3/2 
spins. Therefore, we have five Schmidt coefficients, three of 
which (bosonic part) come from SU(2) 5 = 1 and the remain- 
ing two (fermionic part) come from SU(2) 5 = 1/2. The 
Schmidt coefficients are calculated as 



B. UOSp(l|4) SVBS states 

In the case of UOSp(l |4) ((TV, i^) = (1, 2)), we obtain the 
entanglement spectrum of the MPS (|33|l as: 



Ab^ = Ai^ — X-i' — Aa^ — 



Af^ = ~ — 



5(4 + V25 + 24r2) 



6 6(25 + 24r2 + 4^25 + 24r2) ' 

(57a) 

5(4 + V25 + 24r2) 



4 4(25 + 24r2 + 4\/25T2¥2) ■ 



(57b) 

The bosonic part is triply degenerate as in the case of origi- 
nal 5 = 2 VBS chain, while the fermionic part, which newly 
appeared in SUSY case, is doubly degenerate. Such double 
degeneracy is a fingerprint of a symmetry-protected topologi- 
cal (Haldane) phase in ID.'^ In the absence of fermionic holes 
(r = 0), the fermionic part of the spectrum is infinitely higher- 
lying (see Fig.|6]l and the entanglement of the system is com- 
pletely determined only by the bosonic part which does not 
show the signature of the Haldane phase. 

In the SUSY case, on the other hand, the fermionic lev- 
els appear above the finite entanglement gap and there always 
exists doubly degeneracy in the Schmidt coefficients which 
accounts for the topological stability of the SVBS state re- 
gardless of the parity of the bulk superspin S. We will revisit 
this in section |V] As shown in Fig|6] the five Schmidt coeffi- 
cients take the same value 1/5 at r = 5, and the asymptotic 
behaviors of the entanglement entropy are 



lim 5ee(7') = log 3, 

r— i-O 

lim 5EE(r) =log2+ ;^log6 

r— >-oo z 



(58) 



Thus, at r — > oo, the 5 = 2 SVBS state supports the finite 
entanglement entropy and does not reduce to a simple product 
state as in the 5 = 1 SVBS chain. 



(A.(r))^ 
(A5(r))^ 



1 5 

8 8V'l6r2 + 25 

1 5 



(a = 1,2,3,4) 



(59) 



2^/T6^ 



25 



which are plotted in Fig. |7] together with the correspond- 
ing entanglement entropy. The bosonic part of the spectrum 
is quadratically degenerate while the fermionic part is non- 
degenerate. In both cases, the entanglement entropy 5ee('') 
takes its maximal value log 5 at intermediate value of r = 5/3 
where all the five Schmidt coefficients coincide. The entan- 
glement entropy 5ee('') exhibits the following asymptotic be- 
haviors: 



lim 5ee(?') 



lim 5ee(7') = log 4 



(60) 



If we had a boson instead of the fermion in ( [33] ) as in 
the boson-pair VBS state eq.(|56]l, entanglement would vanish 
in the limit r ^ oo. Therefore, the existence of finite entan- 
glement even in the r oo limit may be attributed to the 
fermionic property of the holes. 

Here it should be emphasized that all the limiting behaviors 
( |54l ), ( |58] l and (|60| can be understood from the viewpoint of 
the edge states; basically, the limiting value of 5ee('') is deter- 
mined solely by information of the irreducible representation 
which describes the emergent edge states. In fact, the general 
formulas ( |B9[ ) and ( |B12| i given in appendix [B] reproduce the 
above results. 



V. SUPERSYMMETRY-PROTECTED TOPOLOGICAL 
ORDER 

In this section, we show that a family of SVBS states 
|SVBS(l|2i^)) exhibits the generalized topological order 
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FIG. 7. (color online) The entanglement spectrum of the U0Sp(l|4) 
SVBS state ([37j. The entanglement entropy of the same state is 
shown in the inset. 



which will be characterized below. Our argument is a SUSY 
generalization of the one presented in Ref. [iTl In the follow- 
ing arguments, we utilize the SMPS formalism. The SMPS 
formalism itself is defined independent of the super Lie group 
symmetries, and is a general formalism to treat a system 
of boson-fermion mixture whose ground state is represented 
by a supermatrix. Therefore, though we mainly discuss the 
SVBS states which have specific underlying particular super 
Lie group symmetries, the following arguments apply to any 
boson-fermion mixture systems. 

Before going into the detail, we first characterize the sym- 
metry operation (both unitary and anti-unitary) within the 
framework of MPS^ The MPS A1A2 • is said to be in- 
variant under the (anti-)unitary operation if the transformed 
state A'iA2 ■ ■ ■ coincides with the original one up to an over- 
all phase. Then, it can be shown"*^ that the invariance of a pure 
MPS is equivalent to the existence of a I?-dimensional {D be- 
ing the size of the MPS matrix A) unitary matrix U which 
satisfies 



(61) 



The phase 9 is not universal and depends, in general, on the 
symmetry operation in question. 

The (c-number) unitary matrix U in ( |6T] l may be postulated 

as: 



Ui 



Up 



(62) 



where Ub and Up are unitary matrices that act on the two 
bosonic subspaces having different fermion numbers. The 
reason for choosing the above form may be seen as follows. 
First we note that eq.(|6T]l implies that the MPS transforms like 



\'f)^Str{UUiA2---A2n+lU) , 

where supertmce is defined as 



str 



A\^) /l(2) 



1(2) 



= trA«-tr4^) 



(63) 



(64) 



While in the case of bosonic MPS, this, combined with 
lx{AB) = tr(B^), immediately implies 11*) cx |*), the rela- 
tion str(Ai3) = str(i3yl) holds only when A and B are super- 
matrices (that contain the Grassmann-odd blocks in their off- 
diagonal parts). In fact, if A and B were merely the c-number 
matrices, A and B, in general, would not commute inside 
str(-): str(AB) ^ str(Byl). To satisfy str(AS) = str(SA) 
only with c-number matrices, either ^ or S is forbidden to 
have c-number components in the off-diagonal blocks. 

Physically, the above relation states that the original sym- 
metry operation (acting on the physical Hilbert space on each 
site) 'fractionalizes' into the ones {U and U^) which act on 
the edge states on both ends of the system. 

In what follows, we parametrize the yl(rn) -matrices in 
terms of the Z) x D-matrices (A, {r(TO)}) as A{m) — 
r(m)A. The diagonal matrix A contains the Schmidt eigen- 
values in its diagonal elements (tr(A^) = 1) and commutes 
with the unitary matrix: [A, U] = 0. In what follows, we use 
the symbol F for the MPS A-matrices in the canonical form.22] 

Then, the F-matrices satisfy the condition for the canonical 
MPS on infinite-size system^ 



^r^(m)A^F(m) = Id 



(65) 



(For more details about the properties of U, see appendix [C]) 
In terms of these F matrices, eq.(|6T]l reads 



F'(m) = e*V^F(m);7 



(66) 



Now let us determine the properties of U satisfying the above 
equation for specific symmetry operations. 



A. Inversion symmetry 

A matrix product state on a circle is given by 

1^-) =str(^i^2---^2„+i) , (67) 

where 'str' denotes the super-trace. By the inversion with re- 
spect to a given link, the state is transformed as 



11*) =Stl-(^2n+l---^2^l). 



(68) 



Here, we use the property of the supertrace: str(MiAf2) = 
str((MiM2)'') = str(M;^'Mf ) to rewrite the above as I|*) 

Z|*) -str(^lUf •••^^'„+i) , (69) 

where supertransposition 'st' is defined as 



(70) 



Therefore, the link-inversion I amounts, in terms of A, to 



(Ml 




( M{ 




\N2 


M2 - 




M'2) 



If we write 



A = ^AF(m)|7 



(71) 



(72) 
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we see that X acts on r(TO) as 



r(m) A r'(m) = r(m)^ 



(73) 



Here, m labels both bosonic and fermionic components and 
r(m) are given by 



r(m) ^ ( ^■^1(5"') ^^^0^^) (m: bosonic) 

Vim)— ( ,r ? \ (to: fermionic) 

\jV2(to) / ^ ' 



(74) 



Originally, Mi , , A^i and N2 are all c-number coefficient 
matrices. However, for practical reasons, it is often convenient 
to assume that the basis states are commuting and take into 
account the anti-commuting properties of the fermionic states 
by supermatrices. 

If I leaves the MPS invariant up to a phase, the general 
relatiorP^ (|66| implies that there exists a unitary matrix [// 
satisfying 



(75) 



In fact, we can prove that d can take the only two values, 
and TT, namely 



L//^r(m)[// = ±r(TO)' 



(76) 



For later convenience, we introduce the following diagonal 
matrix having the same block diagonal structure as C//: 



P 



1b 
-If 



(C//P = PUi) . 



(77) 



Then, the fact that the link-inversion squares to unity leads 
to an important conclusion that Ui is a 'symmetric' or 'anti- 
symmetric' unitary matrix: 



(78) 



The appearance of P is closely related to the property of su- 
pertransposition: 



(A''f = PAP 



(79) 



We give the outline of the proof in the appendix [C] 

By computing the determinant of the above, one can show 
that either fermionic (when the sign + occurs) or bosonic (— ) 
sector has even-fold degeneracy in each entanglement level, 
which we will use as the fingerprint of the SUSY-protected 
topological order. 



case. Under the time reversal transformation T, the spin is 
transformed as 



T 



-Sa- 



(80) 



In the usual matrix representation, the above relation can be 
expressed as 

Sa ^ -Sa = (e"^»X)5a(ife-"^«) = Rl,{7r)St, (81) 

where K is the complex conjugation operator and i?^(7r) rep- 
resents the TT-rotation around the y-axis: 




(82) 



As in the usual case, time reversal operation is defined as 



Sa A {e'^^yK)SaiKe-'^^y) = ~Sa 



S^ ^ (e'^'^y K)S„{Ke 



— i-rrSy 



(83) 



where UOSp(l|2) superspin matrices Sa {a = x, y, z) and S^ 
{(J ~ 61,92) are defined as 



Sa — 



1 



aa 




S^ = 



1 





{ia2r,y 



(84) 



with the Pauli matrices CTq. and ti — (1, 0)' and T2 = (0, 1)'. 
The fermionic generators S^ have the off-diagonal blocks 
which transform as different irreducible representations of 
SU(2) and act as spin-1/2 raising- and lowering matrices. In 
the Schwinger operator representation, S^ are explicitly given 
by Sg, = i(at/ + pbl Sg, = i(6t/ - /ta). Under the 
time-reversal transformation, the SU(2) spinor states are in- 
terchanged: It) = a^|0) ^ \i) = 6t|0), li) = 6^|0) -> -It 
) = — a^|0), and the spin-less fermion state remains the same: 
/^|0) — > /^|0). This implies that the time reversal transfor- 
mation of Sa is given by (83 1. Then we have T^Sa — —Sa, 
so the relation = — 1 for half-integer spins appear for the 
"fermionic spins". 

In fact, for integer superspins, T satisfied 



T^^V, {V)r 



-1) 



F{n) 



(85) 



where V acting on the physical Hilbert space is analogous 
to P in eq.(|77|i acting on the auxiliary space and, due to the 
fermion number operator F{n) {F{n) = or F{n) ~ 1 when 
n labels the bosonic or fermionic variables), (— 1)^(") gives 
a minus sign for the fermionic sector of the (physical) Hilbert 
space. 

Using the above properties, one can readily see that the time 
reversal operation transforms r(TO) as: 



(86) 



B. Time-Reversal Symmetry 

Before discussing the properties of SMPS under time- 
reversal, let us define the time-reversal operation in the SUSY 



Then, time reversal invariance of the SMPS means that there 
exists a unitary Ut such tha^ 

^ i?L.(^)r*(n) = e'^-4r(m)f/T, (87) 
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withe'^^ = ±1. The property = P (for integer superspin) 
requires that the unitary matrix Ut should satisfy 



U't = ±PUt 



(88) 



Since this is exactly the same as eq.(78 i for the link-inversion, 
a similar conclusion is drawn about the entanglement spec- 
trum. 



C. Z2 X Z2 symmetry 

The Z2XZ2 symmetry*^ in the original bosonic case is gen- 
erated by the two commuting tt rotations around x- and z 
axes. However, the symmetry around each axis alone does 
not directly imply the double degeneracy of the entanglement 
spectrum. Rather, it has been shown'' that their combination 
leads to a non-trivial conclusion concerning the entanglement 
spectrum. Below, we show that an analogous symmetry leads 
to a similar conclusion even in the presence of SUSY. 

The TT rotation around the x (z) axis Ux{t^) {Uz{t^)) acts on 
SMPS as: 



r(TO) °^ \ r(m)' 



i?^„(7r)r(n) ia^x,z), (89) 



where _R5^„(7r) is the (45 + 1) -dimensional rotation matrix 
of UOSp(l|2) (see, e.g., eq.( |C^ ). The right hand side is 
equivalent to the action of a unitary matrix C/J^ 



Y.R'LnMnn) = e*'^"C/tr(m)C/, (a = x,z) 



Then, the property = V implies the following 



g2,:e. = 1 ^ 6^*^== = ±1, 



(90) 



(91) 



The phase factor e**"" can be absorbed in the definition of Uj- 
and we may assume Ul = PUa (a = x, z) hereafter. 

On the other hand, for the combination of the rotations 



Ma;(7r) and Uz{tt), we obtain (see appendix C 3 for detail) 



(92) 



e'"-' = ±1, 

By using Ul ~ PUa obtained above, one can show e* 
± 1 and the following exchange property: 



U^Uz = ±PUzU, . (93) 

In terms of the block components Ua,B and Ua.F, this reads: 

C/,,Bf/z,B = ±Uz.bUx,b , Ux.fUz,f = tUz,fUx,p , (94) 

which immediately implies the same degenerate structure of 
the entanglement spectrum as in the two previous cases. 



D. (Z2 X Z2)^ symmetry in UOSp(l|4) SVBS 

Now let us discuss the entanglement spectrum in the 
systems with SO(5)-symmetry and its SUSY generalization 
UOSp(l |4). Inversion symmetry acts independently of the in- 
ternal symmetry and leads to exactly the same conclusion as 
above. The crucial difference from the SU(2) case is the exis- 
tence of (Z2 X Z2)^-symmetry^'* in a class of the SO(5) VBS 
states!^ Specifically, the group (Z2xZ2)^ consists of the fol- 
lowing 16 elements: 



X (l,i?'^(7r)) X (l,i?^*(7r)) X 



(95) 



with R°'^{'k) = exp(i7r(Tafc) iaab'- SO(5) generators). The 
four-fold degeneracy of the entanglement spectra of the SO(5) 
VBS states has been discussed^ from the viewpoint of (Z2 x 
Z2)^-symmetry. 

It is straightforward to generalize the above symmetry to 
the UOSp(l|4) case; now the matrices R"-''{tt) satisfying 
R°'^{'k) = 1 are replaced by the block-diagonal matrices of 
the fomP 



R 



(B) 







R 



(F) 



(96) 



For instance, in t he su perspin-1 UOSp(l|4) SVBS state dis 
cussed in section 



lie 



i?^^' and R^2 ^"^ given by i?'^''(7r) in 
the adjoint- (10) and the spinor (4) representation of SO(5), 
respectively. It is easy to show that the above matrices satisfy 



i?°^(7r)i?"^(7r 
R^\n)R'\TT 



7^4110 (no sum for a and 6) 

7'4|ioi?''(^)i?''(7r) 

R^\ir)R^\^) , 
R^HTT)R^^{Tr) , 



(97a) 
(97b) 
(97c) 

(97d) 



with 



V4 



|io 



lio 

-I4 



(98) 



Now we can apply the argument in section V C 



smce we 

have the same exchange relations ( |97a| l, ( |97b 1 as before. 
Then, we immediately conclude that there exist two sets of the 
corresponding unitary matrices {Ui2,U25} and {1/34,1/4^} 
satisfying 



U42U25 = ±PU25U42 , U34U45 = ±PU45U34 , 



(99) 
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where the matrix P is defined in eq.fTT]). Note that the same 
sign should be chosen for the two exchange relations by the 
SO(5) symmetry. 

The role of the unitary transoformation Uab is clear. First 
we note that, as in the SO(5) case, the following two are mu- 
tually commuting generators of the same block-diagonal form 
as R''''{n) (eq.^) 




(100) 



and can be used as the weight of UOSp(l |4). Since R^^ and 

R'^^ act on the weight (L^^, L^^^) as 



^25t2^12^25 ^ R^^^ L^"^ ~ 



(101) 



uating them using MPS, the following matrices are necessary 



m.n— 1 

d 



m,n— 1 
d 



\rpa 1 _ 

L string] a, q;/3,/3 



(a — X, z) 



mlS'^e"^" 



(103) 



as well as the usual transfer matrix. For instance, the MPS 
expression of the string order parameter Oftring (f^i" open 
chain) reads: 



it is legitimate to assume that the algebra is represented in 
the product space Vi®V2 where Vi and V2 respectively cor- 
respond to {[/12, U25] and {1/34, C/45}. For instance, the two 
unitary operations C/25 and C/45 actually mean 



string 



51 exp 



j+n-l 

in J2 SI 

k=j 



j+n 



^ -'string I-' String J -L -L ; 



(104) 



f/25 «) 1 , 1 «) C/45 ■ 

(C/25<»l)(l®t^45) 
(l®C/45)(t/25® 1) 



U2b ® U45 
(PC/25) ® Ui5 



(102) 



Now we use the fact that Vi and V2 should always have even- 
dimensional sectors and (they have the same dimen- 
sions by the SO(5)-symmetry) to show that the dimension of 
Vi 1S1V2 should be integer-multiple of four This explains 
the existence of the four-fold-degenerate entanglement level 
in the UOSp(l|4) SVBS states (see also the argument in ap- 



pendix C 4 1. 



where we have omitted the denominator necessary to normal- 
ize the MPS. The two parts T^^ (A^l = j - 1) and T^" 
(/Vr = L — n — j) are straightforward; for the canonical MPS, 
they reduces, in the infinite-size limit, to: 



[-^ \a.»:.l3R,PR > Oa^aOp^^fj^ 



(105) 



The boundary dependent factors (^ql.ql and Sj^^ are can- 
celled by those coming from the denominator. Therefore, all 
we have to compute is the infinite-distance limit (n /• 00) of 
the following quantity: 



VI. RELATIONS BETWEEN STRING ORDER 
PARAMETER AND TOPOLOGICAL ORDER 



^ ^ [^String (^string ) T ]a,a;^,/3 
a,/3 



(106) 



Later, the use of the string order parameters in detecting 
the Haldane phase was criticizecf'^ since they are well-defined 
only in a restricted class of models and fail to capture the ro- 
bustness of the Haldane phase as a symmetry-protected topo- 
logical phase (see Refs.|48jand|49jfor the attempts at alterna- 
tive order parameters). Now a natural question arises; under 



what conditions the string order parameters (45ai and (45b 1 



correctly capture the topological nature of the Haldane phase? 
Below we will uncover the explicit relationship between the 
string order and the topological order to answer to this ques- 
tion. 



A. String Order Parameters in MPS Framework 

Let us fi rst consider the structure of the string order param- 
eters (45a I and (45b 1 from the MPS point of view.^"*^ In eval- 



B. String Order Parameters and Entanglement Spectrum 

Now we show that the existence of non-vanishing string 
order parameters serves as the sufficient condition for the 
symmetry-protected topological order discussed in the previ- 
ous section. Let us begin with the simpler case of the usual 
VBS states. 

Since we are interested in the long-distance limit |i — j| Z' 
00, we need to know the asymptotic behavior of the string 
(Tstring)'^"^'- To this end, we can borrow the results of Ref.|43] 
(Theorem 2); according to the theorem, the MPS should be 
invariant under both of the 7r-rotations 



,-i7rS" 



(107) 



in order for the string (Tsning)'* not to vanish in the long- 
distance limit. Then, Lemma 1 of Ref. 43 guarantees that there 
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exists a pair of unitary matrices Ux and Uz which are unique 
and satisfy: 



E R[p{nUnA{n) = e*«" UlA{m)U, 

n=l 

(a = a;,z; e*^° = ±1) 
{Uaf = 1 , UxU, = ±UzUx , 



(108) 



where the two sign choices are independent. The above ex- 
change property between and has a very important im- 
plication to the structure of the entanglement spectrum'^; 

det {(C/,C/,)a} - det {(C/,)a} det {(;7,)a} 

= (il)"^- det {{Ux)x} det {(f/,)^} 0) . 

(109) 

Therefore, the degree of degeneracy d\ of each entanglement 
level A should be even when Ux and Uz are anti-commuting. 
Typically, this happens in the VBS states with ot/cZ-integer-S'. 

Now we show that when the string order parameters are 
non-vanishing O^^^ ^ 0, the minus sign realizes (i.e. Ux 



and Uz anti-commute) in eq.( 109 1 and the entanglement spec- 
trum has the degenerate structure. To this end, we investigate 
eq.( 106 1. First of all, the invariance of the MPS under u^. ^ 
implies that the string part (Tstri„g)"^^ reduces essentially to 
a phase (e'®")"^^ = (±1)"^^. This is a direct consequence 
of Theorem 2 of of Ref. 43 and is easily understood since the 
overlap (^'|-Ua|\l') = (Tltring)^ vanishes otherwise. The price 
to pay is the boundary factors appearing at the two end points 
of the string correlation functions (see FigjSli: 



{ ^string I ^'^Ln'^Ln | V-^ string, 
a,P y \n=l 

|j-j|/'oo 



E{(^-.ngVi:j)(vSr^)} 

^ — ' 1^ ) a, a 



a,/3 



where V^^j^ denotes the left (L) and the right (R) eigenvec- 
tors of Tstring. 

To see whether the boundary factors are non-vanishing or 
not, we consider the right-boundary factor (1 {l®Uz} T^) of 
O J^iij^ (i.e. a — z). First we rewrite it by using (see the second 
figure of Fig|9]l: 

S'' = uluxS'-ulux = ul{-S')ux {ux = «)fee-"^°') . 

(Ill) 

The unitary operators u\. and Ux appearing on both sides of 
— iS*^ can be absorbed into the MPS matrices by using eq.( 108 i 
(the third figure of Fig|9]l. By re-arranging the unitary matri- 
ces UxUz (the fourth figure of Fig|9]l, we arrive at the expres- 




(110) 



l{l®Uz}T'- - 1 {m{UxUzUl)] (-T^) 

= i{m(±UzUxUl)]{-T') 

= T'^{mUz}T' . 



(112) 



Therefore, we see that the boundary factors, and hence the 
string order parameter itself, vanish when Ux and Uz are com- 
muting (like e.g. in the even-5 VBS states). On the other 
hand, if both of the string order parameters are finite, this im- 
mediately implies that the ground state MPS is not only in- 
variant under the two 7r-rotations"^^ Ux and Uz but also has the 
adjoint Ux,z matrices satisfying 



UxUz = -UzUx 



(113) 



By the argument in Ref. IJ the ground state is topologi- 
cally non-trivial in the sense that each entanglement level is 
even-fold degenerate. Therefore, the finiteness of the pair of 
string order parameters O^tVlng 1^ '■1'^ sufficient condition for 
the topological phase. It is crucial that both O^rfng ^string 
are non-zero for the existence of the topological order. For 
instance, one can construct a solvable spin-1 model^'^ which 
exhibits a kind of 'hidden order' similar to the one in the VBS 
model and has^' *^Srinc ~ and 0^^^^^ 7^ 0. In fact, in this 
case, the two entanglement eigenvalues are no longer degen- 
erate and the state is not topological. 




FIG. 8. (color online) Diagrammatic representation of the main part 
of string correlation function | (T,fn„g V^"' ) ( V["i' ) | . 







J; ^ 











-s 
u, u* 



D 




-0 



-D 



FIG. 9. (color online) Rewriting the boundary factor (for a — z) 
using Ux- When Ux and Uz are anti-commuting, the minus sign 
coming from UxS^ul: = —S^ is cancelled and an overall plus sign 
is recovered. 
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C. CaseofSMPS 



VII. SUMMARY AND DISCUSSIONS 



Basically, we follow the same line of arguments to show 
that finite string correlation implies the topological phase. The 
only difference is that now we have the P matrix ( |77| ) in the 
key equation ( 113 i: 



(114) 



Correspondingly, the last step (see Fig. [9]l in evaluating the 
boundary factor is modified. Specifically, in stead of eq.( |l 12} , 
we have (see Fig.[TO|: 



(115) 



i{mu,}T' ^i{m{u,u,ul)}i-T'-) 

= T'i-{mPU,}T' . 



Therefore, one of the two components (bosonic and 
fermionic) vanishes just by symmetry: 




(116) 



when e^'^''' = +1 



when e*'^="^ = — 1 



Therefore, if the two string order parameters are both non- 
vanishing, either the bosonic- or the fermionic sector ex- 
hibits the degenerate structure mentioned in section [V] and the 
ground state is topologically non-trivial. 

Now it is straightforward to generalize the above argument 
to the case of UOSp(l |4) to show that when all the four string 
order parameters 



ab 



\i-]\y^oc \ 



exp 



i-i 



k= 



- ab 



(117) 



(where (a, b) ^ (1, 2), (2, 5), (3, 4) and (4, 5), and Lab are the 
SO(5) generators) are non-zero, 2^ x (integer)-fold degeneracy 
occurs in some (bosonic or fermionic) sectors of the entangle- 
ment spectrum. 



-O 



-s 



D 




-D 



-D 



FIG. 10. (color online) Rewriting the boundary factor (for a = z) 
using u^. When Ux and Uz are anti-commuting, the minus sign 
coming from UxS^ul. = — is cancelled and an overall plus sign 
is recovered. Note that an extra P-matrix appears in SUSY case. 



We investigated the effects of doped fermionic holes on the 
topological phases in quantum antiferromagnets. To this end, 
we first introduced a family of SVBS states which may be 
thought of as the hole-doped version of the usual (bosonic) 
VBS states e.g. spin-S" SU(2)- states, the SO(5)- and the 
Sp(iV) VBS states. One of the standard ways of looking at the 
topological properties in these states is to investigate the string 
order parameters. We explicitly evaluated the behaviors of the 
string order parameters of the UOSp(l |2)- and the UOSp(l |4) 
SVBS states for various values of superspin-5, and found that 
even when the string order parameters vanish identically in 
the absence of doping, they revive immediately after holes are 
introduced in the system. This might suggest that the doped 
holes changes the property of the ground state and thereby 
stabilizes the topological phase. 

To better understand the nature of the states, we calculated 
the entanglement spectrum. Basically, the spectrum consists 
of the bosonic sector and the fermionic one; at zero doping 
r = 0, the fermionic sector is separated from the bosonic 
sector, which constitutes the low-'energy' part of the spec- 
trum, by an infinitely large entanglement gap. Upon doping, 
the fermionic sector starts participating in the entanglement. 
The point is that the existence of supersymmetry allows the 
coexistence of the two sectors having different entanglement 
structures. In addition to that, the entanglement spectra in the 
SUSY systems exhibit the following peculiar features: (i) In 
contrast to naive expectation, the SUSY entanglement spectra 
for the bosonic- and the fermionic sectors do not coincide with 
each other at r = 1, as a consequence of SUSY many-body 
effect, (ii) in the two extreme limits of the doping parameter, 
r and oo, the entanglement spectra of the SVBS states 
indeed reproduce those of the original bosonic VBS state and 
the Majumdar-Ghosh-like states, respectively. 

On the basis of the observations made for the particular 
states (UOSp(l|2) SVBS and UOSp(l|4) SVBS), we charac- 
terized, with the help of the SMPS formalism, the symmetry- 
protected topological orders in the SUSY systems in terms of 
the entanglement spectrum. According to the results, there al- 
ways exists a topologically-protected sector (whose degener- 
ate structure depends on the symmetry of the SMPS in ques- 
tion) in the spectrum of the SUSY systems. Also, by using 
the SMPS formalism, we clarified an intimate connection be- 
tween the finiteness of the string order parameters and the de- 
generate structure of the entanglement spectra; the finite string 
order is the sufficient condition for the degeneracy in the en- 
tanglement spectrum, which is the fingerprint of the (topolog- 
ical) Haldane state in the bulk. These explain the revival of 
the string order upon doping. 

The above remarkable features can be understood in the 
light of SUSY edge state picture. Intuitively, the degener- 
ate structure can be understood by the existence of fictitious 
'edge' superspins that appear at the entanglement cut of the 
chain. When the bulk system has superspin S, two super- 
spins iS/2s, which consist of the SU(2) spin 5/2 and its super- 
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partner 5/2 — 1/2, emerge at the edges: 



ACKNOWLEDGEMENT 



S/2 



SUSY 



S/2 - 1/2 



(118) 



Then, there always exist half-odd-integer spins at the edges 
regardless of the parity of the bulk superspin, since SUSY, be- 
ing the symmetry that relates the state with integer spin and 
that with half-odd-integer spin, guarantees the coexistence of 
both. Such half-odd-integer 'edge' spins bring the even-fold 
degeneracy to the entanglement spectrum of the UOSp(l|2)- 
symmetric systems. Therefore, if we have a topological phase 
(e.g. Haldane phase) characterized by the above type of de- 
generate structures in the entanglement spectrum, it exists for 
all values of superspin S. A similar argument applies, with 
due modification, to cases with other types of SUSY. In this 
sense, one may say that SUSY plays a unique role in stabiliz- 
ing the topological phases of matter in ID. 

Since our study presented here is restricted to a particular 
class of VBS states with SUSY, one obvious future direction 
would be to extend it to more generic models. The argument 
for symmetry-protected topological orders presented in this 
paper can be generally applied to any system whose ground- 
state wavefunction is given by the (S)MPS states. Thus, it 
would be interesting to see, for instance, the robustness of the 
Haldane phase in the SUSY Heisenberg model with respect to 
the parity of the bulk superspin S. This might highlight the 
unique behavior of SUSY topological phases in comparison 
to the bosonic counterparts studied in Ref.llTl 

Another future direction is the generalizations to higher di- 
mensions. In higher dimensions, the SVBS states generally 
interpolate between the bosonic VBS states and the RVB type 
of states (in ID, we have the Majumdar-Ghosh valence-bond 
crystals). The latter is well-known to have non-trivial topolog- 
ical properties and it would be interesting to study the change 
in the entanglement properties and the edge-state structure as 
the doping is varied. Application to other topologically non- 
trivial states of matter, such as quantum Hall states or various 
topological states in cold atom systems, is even more interest- 
ing. For instance, the SUSY-extended Laughlin wave func- 
tion, which has a close analogy with the SVBS states stud- 
ied here, interpolate between different quantum-Hall ground 
states, such as the Laughlin states and the Moore-Read Pfaf- 
fian states. In this respect, as the SVBS states in ID provided 
a unifying way of deriving the entanglement spectra of the 
(bosonic) VBS state and the MG dimer state, the study of 
the entanglement spectra of the SUSY Laughlin wavefunction 
will naturally give a unifying understanding of the entangle- 
ment structure of various quantum Hall ground states. 

Last, we would like to comment on the recent work on the 
non-local order parameters for the symmetry-protected topo- 
logical order. When completing this paper, we became aware 
of a recent preprint by Pollmann and Turner (Ref. 49) which 
also discusses the string order parameter from the entangle- 
ment point of view. Although some of the conclusions ob- 
tained there overlap with ours, the main goal there is to go 
beyond the string order parameter and is different from that of 
the present paper. 
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Appendix A: A-matrices for UOSp(l|4) SVBS states 

1. Superspin-1 SVBS 

The fourteen 5x5 A-matrices for the 5 = 1 SVBS state 
discussed in section [iTC] are expUcitly given as: 

0^ 

= -A(2,2)' = -\/2 O2 
\ Oy 

/O2 0^ 

yl(3, 3) = -A(4, 4)' = -V2 cr_ 

\ Oy 

^0-3 o^ 

A(l,2) = I O2 
0; 



A(l,3) = -A(2,4)' = - 



A(l,4)=^(2,3)' = 



^02 0> 
A(3,4) = ( as 
0; 



02 








O2 












O2 




-12^ 


i + 




O2 











(Ala) 



A{1) = A{2f = 
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(Alb) 



where the symbols 't' and 'st' denote the transposition and 
supertransposition ( 70 1, respectively. They can be represented 
by linear combinations of the UOSp(l |4) generators. 
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2. Properties 



As has been discussed in section |V A| the link-inversion 
symmetry is implemented in the SMPS as 



I : Aim) ^ A{my' , 



(A2) 



or to write the bosonic- and the fermionic component sepa- 
rately 

I : A{a, t) ^ A{a, t)' , yl(cr) ^ A{af . (A3) 
Then, it can be shown 

A{mf ^ yV^A(m)W , (A4) 

where 

'W 0^ 



1 



with 



W 



( ^f72 
\i(J2 



(A5) 



(A6) 



Appendix B: Edge States and General Asymptotic Beliavior of 
Entanglement 

The asymptotic behaviors eqs.(|54]i, ( |58| ) and (|60| can be 
understood from a more general point of view. Let us con- 
sider the UOSp(l|2iir) SVBS state with bulk-superspin S. 
The UOSp(l|2A') SVBS has = 1 supersymmetry, and con- 
sists of one bosonic sector and one fermionic sector For the 
bulk-superspin S, the emergent superspin-5/2 objects appear 
at the edges and the UOSp(l|2i^) SVBS state accommodates 
the graded fully symmetric representation^ at each edge: 



|TOi,m2, • • • ,m,2K) 



(6lr^(4)"^...(5^^r-|vac), 



(Bla) 



|n-i,n-2,--- ,"-2k) 



Vni!n2! • • • ^2^! 



(Bib) 



with mi + m2 + • • • + ■m2K — ni+n2 + - ■ ■+ n2K + 1 = 5. 
Then, the number of the bosonic- and fermionic states on each 
edge are respectively given by 



S + 2K-l\ _{S + 2K- 1) 



D 



' S ) {2K-l)\S\ ' 
5 + 2/^-2\ {S + 2K-2)\ 



S-1 



{2K - 1)\{S - l)V 



(B2a) 
(B2b) 



(The bosonic degrees of freedom coincide with the fully sym- 
metric representation of USp(2A')— .) For instance, for the 
UOSp(l|2) {K = 1) SVBS state, we have 



Db = 5 + 1 , Dp = 5 



(B3) 



while for the UOSp(l |4) {K = 2) SVBS state, 

= ^(5 + l)(5 + 2)(5 + 3), 

D 

D^=\s{S + l){S + 2). 
6 



(B4) 



In the infinite chain limit, the spin degrees of freedom are 
equivalent 



Ai^ — _ . . . _ \db^ — Ab^, 
A_Db+i^ = Adb+2^ = • • • = a 



(B5) 



and the normalization condition of the Schmidt coefficients, 
Sa=-i^'' I'^aP = 1' is rewritten as 



Db-Ab^ + i?F-AF^ = 1. 



Then, the entanglement entropy is expressed as 



(B6) 



SeeM = - ^ |A„plog|A„p- ^ |Az,B+aplog|Az,B+„|2 

a— 1 a—1 
= -i?B|ABplog|ABp-i?F|AFplog|AFp. 

(B7) 

At r = 0, only the Schmidt coefficients of boson sector sur- 
vive and eq.(|B6[) implies 



and hence 



Ab^ — 7:;-, Af^ — 0, 



lim S'ee('') = log 15b • 



(B8) 



(B9) 



Thus, the entanglement entropy of the spin 5* original VBS 
states is reproduced. 

On the other hand, in the limit r ^ 00, the SVBS states 
reduce to the (partially) dimerized states [see FigfTT]. 



+ 



o o 



4^ 
4}- 



FIG. 11. (color online) The r —> 00 limit of the bulk superspin <S = 
2 U0Sp(l|4) SVBS state is given by the superposition of the two 
partially dimerized states related to each other by 1-site translation. 
When we make an entanglement cut at an arbitrary bond (shown by 
wavy lines), we always have two different kinds of sections: one with 
four 'fermionic' edge states (upper) and the one with ten 'bosonic' 
edge states (lower). These two different sections respectively yield 
four-fold- and ten-fold degenerate entanglement levels. 



In the upper state in Fig. |2j the fermionic edge states ap- 
pear, while in the lower the edge states are bosonic. Since 
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both cases appear with equal weights, the sum of the Schmidt 
coefficients for the bosonic sector and that for the fermionic 
sector should be equal: 



-Df 



Aq^ — ^ Xdb+o^ — 1/2. 



Therefore, we have 



A 



B 



1 

2np 



(BIO) 



(Blla) 



for r — > oo, and the corresponding entanglement entropy is 
derived as 

lim 5e.e.('^) = log (2V^b^f) , (B12) 



with I?B and Df given by eq.(B4i. Thus, from the entangle- 
ment point of view, the role of SUSY is two-fold. First, it 
necessitates two different Schmidt eigenvalues corresponding 
to the iV = 1 SUSY. Second, it enables the system to support 
finite entanglement even in the limit ?- — > oo. 

For the superspin-5 UOSp(l|2) SVBS states,'23 the entan- 
glement entropy behaves as 



(B13a) 
(B13b) 



lim^EEW = log(5 + l), 
lim 5EE(r) = log2 + i log(5(5 + 1)) , 

r— f oo Z 



which, for 5 = 1 and 5 = 2, reproduces the results ( [54) l 
and (Bsll. For the superspin-5 UOSp(l|4) SVBS states, on the 



other hand. 



limS'EE(r) = log(5 + 1){S + 2)(5 + 3) - log 6, (B14a) 

r— 



lim SEE{r) 



- log3 + log(5 + + 2) + - log 5(5 + 3). 



Setting 5 = 1, we reproduce the previous result ([60] 



(B14b) 



Appendix C: Proofs 

In this appendix, we outline the proof of the important re- 
lations ( |78] l, ( |88] l and ( |93] l. For later convenience, we derive a 
useful property of pure canonical MPSs. 

Suppose that we have a pure MPS whose canonical form is 
chai-acterized by the MPS data^- ^ (A, V) and that it satisfies 
the following relation for some unitary matrix U: 

r(m) = e'^^ U^T{m)U . (CI) 

Since the MPS is canonical, the following holds: 

^V\m)K^T{m) = Id ■ (C2) 



Physically, it states that the £)^-dimensional vector v[°'' 



{Y'l^)a;b = 5ab {l<a,b<D) (C3) 
is the dominant left-eigenvector of the left transfer matrix 

iTLha-rb,b ^ 5](Ar*(m))g5(Ar(m)),fc . (C4) 



Plugging T\m) = e"*""^ Wr^m)U into (|C2]), we obtain: 

^-^0u [/trt(m)C/A2r(m) = Id , (C5) 

m 

or equivalently 

(m) A[/Ar(m) = e^'^" U . (C6) 
This implies that the unitary matrix 

a a 

when viewed as a -dimensional vector, is the left- 
eigenvector of Tl with the eigenvalue e*^" : 



(C7b) 



Since, by assumption of canonical MPS, Id is the unique left- 
eigenvector with the eigenvalue |A| = 1, we conclude 



(C8) 



Since in deriving the above, we have only assumed that the 
(infinite-system) MPS in question is pure and takes the canon- 



ical form, ( |C8) holds for any MPS (including SMPS) satisfy- 
ing the assumption. 



1. Inversion-symmetry 

We use the property = 1 to derive the important property 
(78 1 of the adjoint Ui matrix. Applying supertransposition st 
on (75 1 and using (A'*')'*' = PAP, we obtain 



r(m) 



{UiPUr)^T{m) {UiPUr, 



Postulate U is the block diagonal matrix 



U 



C/f 



(C9) 



(CIO) 



By eqs.(Cl i and (C8i, this implies that the DxD matrix 
{UiPU'j) should be equal (up to an overall phase) to the unit 
matrix: 



{UiPU*i) 



(Cll) 



After multiplying U\ from the right and making transposition, 
we deduce 



Ui 



P^Ul 



Ui <S4> e 



±1 (C12) 



Therefore, we obtain eq.(78 1: 

U} = ±PUi . 



(C13) 



It is interesting to calculate Uj for superspin-S* UOSp(l|2) 
SVBS states. For the 5 = 1 SVBS state, U is identified as 



1 
[// =7^ = I -1 
^0 0-1, 

which satisfies 

ujr{m)Ui = -r(TO)^', 

U} = -PUi . 
For the 5 = 2 SVBS state, U is identified as 



(C14) 



Ur = 



/O 10 \ 

'O -1 ' 

1 
0-1 

Vo 1 y 



(C15) 



(C16) 



and r(m) satisfy the relation 

ujr{m)Ui = +r{my\ 

U\ = PUi. 



(C17a) 
(C17b) 



For the 5 = 1 UOSp(l |4) SVBS state, we use the relations 
given in appendix] A l|to show that 



(C18) 



r(TO)^' = u]v{m)Ui 

U\ = -PUi 



with Ui = W defined in eq.( |A5| l. This is consistent with 
the existence of the four-fold degenerate entanglement level 
in this state (see Fig.|7]i. 



2. Time-reversal symmetry 

If the MPS is invariant under time-reversal, the F-matrices 
satisf>HI] 

53i?^„Wr*(n) = e^^-4r(m)f/T , (C19) 

n 

where the rotation matrix R1^^{tt) takes the block-diagonal 
form 

with Rg{TT) and R^_-^^^^{Tr) being the ordinary rotation ma- 
trices for spin-5 and {S — 1/2), respectively. Since = V, 

(-i)^(') r(/) = X: RL |e R'mnr*{n)\ 

m=l l^n— 1 ) 



m— 1 

{UtU:^}"^ r{i) {UtU:^} , 



(C21) 



or equivalently 

r(0 = {UtU^}^ (-i)^(') r(0 {c/tC/^} 

By using the property 

(-i)^(')r(/) = pr(/)p, 



eq.(C22 1 may be rewritten as: 

r(0 = {c/tPC/t}^ r(0 {UtPU^} . 

Now we can apply eqs.( |Cl| l and ( |C8[ ) to conclude 

U't = ±PUt . 

For 5 = 1 UOSp(l|2) SVBS state, with r(l) 
r(2) = AiO), r(3) = A{-1), r(4) = A{l/2) 
A{-l/2) (|22]i, and 



Ry{TT) 




/O 1 0\ 

'O -1 * 

1 

1 

Vo 



we have 



and 



-1 0/ 

5^i?5^„(7r)r*(n) = 4r(m)t/r 



-PUt . 

For 5 = 2 UOSp(l|2) SVBS state, with 



Ut = 



Ry{TT) 








1 




















-1 





















1 




































-1 












Vo 








1 


0^ 

























1 
























-1 

























1 

























-1 

























1 




















































-1 

























1 

























-1 










\o 














1 








0/ 



we have 



and 



^i??^„(7r)r*(n) = +4rMC/T 



Ul = +PUt 
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3. Z2 X Z2 -symmetry 
Finally consider the tt rotation around the x- and the z-axis, 
r(m) ^ J2 K™('r)r(n) (a ^x,z). (C32) 

n 

Instead of (i?")^ = 1 in the bosonic case, R°- in the SUSY 
case satisfies (i?")^ = T-". Therefore, the use of the terminol- 
ogy 'Z2 xZ2-symmetry' is not precise. However, to underline 
the connection to its bosonic counterpart we use the terminol- 
ogy in the SUSY cases as well. 

When the MPS has such a symmetry, we havd^S 

^<„„Wr(n)=e^«"t/tr(m)t/a (a = z) (C33) 



for some block diagonal unitary matrix: 



_ (Ua,B 

Ua.F 



(C34) 



Now let us consider what (C33i implies. We begin by 
(i?")'^ — V (valid for integer superspin-S): 



{V)nnT{n) = PT{n)P 

= e^''"^i?;^„(7r)C/tr(m)C/, 



(C35) 



{UlfT{n)Ul 



which, after Ps are rearranged, reads 

V{n) = e''''^{UaPUa)^T{n){UaPUa) 

implying 



(C36) 



(C37) 



The phase e*"^" can be absorbed in the definition of Ua and we 
have: 



[UaPUa) ^ Ul = PUa 



(C38) 



Next, we consider the product of the two rotations and 
i?^. In the case of SUSY, they obey the following exchange 
relation: 



(C39) 



When combined with eq.( |C33[ ), this translates into the follow- 
ing relation for F: 

{U^U,yT{m){a,U,) = {U,PU^)^T{m){U,PU^) . 

(C40) 

After rearranging the C/s, we arrive at the form to which 



eqs.lCl I and ( |C8| l are applicable: 

F(m) = {U,PU^Ulul)^T{m){U,PU^Ulul) . (C41) 
Therefore we have 

U,PU^^Ulul^e^"lD (C42) 



with e"^="^ = ±1. The resulting equation 

U^U, = ±PU,a, (C43) 

or 

U,aU,,i = ±C/.,i(7,,i , C/,,2f/.,2 = TC/.,2C/x,2 (C44) 

implies the degenerate structure of the entanglement spec- 
trum. 

Let us calculate [/-matrices for superspin-iS UOSp(l|2) 
SVBS states. For odd-5, they assume the following form: 



^(5-l)/2W 



which satisfy 



(C45a) 
(C45b) 



Therefore, the degenerate spectrum appears in the bosonic 
sector 

For even-5, on the other hand, they are given by: 



t;(5)^^(5/2)(^)^(^^5/^2 



satisfying 



(C46a) 



(C46b) 



which implies that the fermionic spectrum exhibits the degen- 
erate structure. 



4. (Z2 X Z2 ) symmetry 

In this appendix, we summarize some useful relations con- 
cerning the v4-matrices of the UOSp(l|4) 5 = 1 SVBS states 



given in appendix A 1 



The invariance of the MPS under R"-'>{Tr) defined in eq.(|96 
implies'*— the existence of the 5x5 unitary matrices Uab satis 
fying 
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Specifically, Uab ai"e given by 



(C47) 



U12 = 


























— i 






1 



















i 














-1 













— i 




















1 







i 














Vo 















\o 














(C48a) 



U34 = 


















1 

















1 

















-1 





Vo 














C/45 = 








—i 
















—i 





i 

















i 











Vo 














(C48b) 
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It is easy to check that these matrices satisfy 



where 



I4 

-1 



(C50) 



iUi2f = ((725)' = = (C/45)' - ^1|4 

U12U25 = —'Pi\/JJ25Ul2 , — — 7^114 C/45J734 , 

U25U45 = — 7'i|4?745C/25 

U12U34 = U34U12 , U12U45 = U45U12 , U25U34 = U34U25 , 

(€49) — 



By the general argument in section [VP one concludes that in 
some sectors all the entanglement levels are fourx(integer)- 
fold degenerate as is seen in Fig. |7] 
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